Spin dynamics in the low-dimensional magnet TiOCl 
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We present detailed ESR investigations on single crystals of the low-dimensional quantum magnet 
TiOCl. The anisotropy of the ^-factor indicates a stable orbital configuration below room temper- 
ature, and allows to estimate the energy of the first excited state as 0.3(1) eV ruling out a possible 
degeneracy of the orbital ground state. Moreover, we discuss the possible spin relaxation mecha- 
nisms in TiOCl and analyze the angular and temperature dependence of the linewidth up to 250 K 
in terms of anisotropic exchange interactions. Towards higher temperatures an exponential increase 
of the linewidth is observed, indicating an additional relaxation mecharusm. 
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I. INTRODUCTION 

Titanium-based oxides have attracted considerable in- 
terest due to the subtle interplay of spin, orbital, charge, 
and lattice degrees of freedom in these compounds. Es- 
pecially, the tendency of the t2g triplet to orbital degen- 
eracy in nearly ideal octahedral symmetry has triggered 
the search for exotic behavior of the orbital degrees of 
freedom, such as an orbital- liquid state in LaTi03fiiS*2i4 
or the presence of strong orbital fluctuations in YTiOsj^ 

Recently, the fascinating system TiOCl came into fo- 
cus as a possible candidate where orbital ordering in- 
duces a quasi-one-dimensional magnetic behavior and a 
spin-Peierls-like transition to a non-magnetic state below 
Tci = 67 K& The true nature of this transition and the 
magnetic dimensionality of this compound, however, is 
still under debate especially because it is preceded by 
another phase transition at Tc2 ~ 90 K: 
Given the crystal structure of TiOCl which consists of 
Ti-0 bilayers within the (a6)-plane, well separated by 
CI ions (see Fig. QJ,^ Seidel et al. suggested two possi- 
ble chain-directions of the Ti^+ ions, the first one medi- 
ated by superexchange coupling along the a-axis (zigzag 
chain) and the second one by direct exchange along the 
&-axis (linear chain) .5 By LDA-hU and LDA-FDMFT cal- 
culations it was concluded that the direct exchange along 
b dominates the magnetic behavior in agreement with the 
description of the high-temperature magnetic suscepti- 
biliy in terms of an antiferromagnetic S = 1/2 spin chain 
with strong Heisenberg exchange J « 660 K.^-^'^-^^ How- 
ever, as already anticipated by Seidel et al., the observed 
two successive phase transitions cannot be attributed 
to a conventional spin-Peierls transition only, but may 
have to involve interchain couplings and frustration sce- 
narios, which also have to account for the first-order 



characteriiSii^ of the transition at Td = 67 K. Further pe- 
culiar features of TiOCl above Td have been reported by 
NMR, suggesting the presence of several inequivalent Ti 
and CI sites and an incommensurable orbital orderingsS 
Moreover, it was pointed out by Shaz et al. that the sym- 
metry in the temperature range T^i < T < Tc2 is lower 
than the orthorhombic structure at room temperature^ 
Recently, the appearance of these two phase transitions 
has been described within a spin-Peierls scenario as a 
result of frustrated interactions arising due to the lay- 
ered structure of TiOCl»i^*i^ Additionally, the existence 
of strong orbital fluctuations up to 135 K or even room 
temperature has been evoked both theoretically and 
experimentallyjiSiiSiiSiii Previously, electron spin reso- 
nance (ESR) data suggested that there might be signif- 
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FIG. 1: (color online) Crystal structure of TiOCl at room 
temperature (space group Pnma). 
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icant changes in the spHttings of the d-orbitals between 
Tc2 and room temperature, based on the temperature de- 
pendence of the anisotropic g-va\uesm^ In contrast, sig- 
nificant orbital fluctuations have been discarded by op- 
tical spectroscopy studiesi^i and polarization dependent 
ARPES measurementsii^ However, the first excited d- 
level could not be detected by the optical measurements 
leaving the question of a possible degeneracy of the lowest 
lying c?-orbitals unsolved. Using ESR we reinvestigated 
TiOCl in detail and find almost temperature-independent 
g-values up to room temperature in agreement with the 
optical and ARPES studies. Moreover, we discuss possi- 
ble spin relaxation processes in this compound and ana- 
lyze the temperature and angular dependence of the ESR 
linewidth in terms of the symmetric anisotropic exchange 
interaction and the antisymmetric Dzyaloshinsky-Moriya 
interaction. 
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FIG. 2: (color online) Temperature evolution of the ESR spec- 
^^^J^J:'^ ^^/t^^^}!^^^ ^^rum in TiOCl for H \\ a. Solid lines represent fits using a 

Lorentzian line shape. 



EXPERIMENTAL DETAILS 



Single crystals of TiOCl were prepared by chemical va- 
por transport from TiCla and Ti02 ^ The samples have 
been characterized using x-ray diffraction, specific heat, 
and magnetization measurements. The crystal struc- 
ture at room temperature was found to be orthorhom- 
bic (space group Pnma) with lattice parameters of a 
= 0.379 nm, b = 0.338 nm and c = 0.803 nm. The 
magnetic properties were found to be in excellent agree- 
ment with published resultsAi^ The good quality of the 
crystals has been clearly confirmed by the susceptibility 
measurements which reveal a hysteresis at the first-order 
transition at Tci — 67 K which had not been reported 
previously. 

The ESR experiments have been carried out with a 
Bruker ELEXSYS E500 CW-spectrometer at X-band fre- 
quency {ly ~ 9.4 GHz) in the temperature range between 
4.2 and 500 K with continuous gas-flow cryostats for He 
(Oxford Instruments) and N2 (Bruker). ESR detects the 
power P absorbed by the sample from the transverse 
magnetic microwave field as a function of the static mag- 
netic field H. The signal-to-noise ratio of the spectra is 
improved by recording the derivative dP/ dH using lock- 
in technique with field modulation. Dielectric measure- 
ments were performed at temperatures 300 < T < 500 K 
over a frequency range I Hz < < 1.08 MHz using a 
Novocontrol a-analyzer. 



III. EXPERIMENTAL RESULTS 

ESR spectra obtained for TiOCl in the paramagnetic 
regime at different temperatures are displayed in Fig. |21 
The spectra consist of a broad, exchange-narrowed reso- 
nance line, which is well fitted by a single Lorentzian line 
shape. The intensity of the ESR signal /esr is propor- 
tional to the static susceptibility^^ and, hence, exhibits 



also the sharp drop at Td to the non-magnetic state and 
the kink at Tc2- For our samples, Iesr{T) is in agree- 
ment with the results obtained in Ref. 18 and, therefore, 
not shown here. 

The temperature dependent ESR linewidth A_ff and 
the effective g-factor are depicted in Figs. EI a) and|3|[b), 
respectively. As reported previously, both quantities 
show an anisotropic behavior with the external mag- 
netic field H applied along the three crystallographic 
axes of the orthorhombic structure>i^ Above T ~ 90 K 
the linewidth is largest when H\\a while the values for 
H\\h and H\\c are almost equal. The linewidth increases 
monotonously for all three directions for T > Tc2 K, how- 
ever, a peculiar change from a negative to positive curva- 
ture is observed at about 250 K. Below T ~ 90 K there is 
a crossover of the linewidth data resulting in the broadest 
spectra for H\\b, which is highlighted in Fig. [T] On ap- 
proaching the first-order transition at T^i, the linewidth 
for all directions drops down to a value of about 50 Oe. 
This corresponds to the residual signal due to paramag- 
netic impurities which will not be further discussed. Fo- 
cusing on the high-temperature behavior above 250 K we 
find that the anisotropy of the linewidth becomes smaller 
and vanishes at about 430 K. 

Notably, at the same temperature the anisotropy of 
the effective g-factor vanishes, too (see Fig.OIb)), and we 
obtain g(430 K) ~ 1.96 for all three directions. Concomi- 
tantly with the change of curvature of the linewidth at 
250 K the temperature dependencies of the g-factor show 
a steep increase above 250 K, while the g-tensor is nearly 
constant in the temperature range Tc2 < T < 250 K. This 
behavior differs from previously published results where 
a much larger and temperature dependent anisotropy of 
the g-factor was reported for Tc2 < T < 300 K and in- 
terpreted in terms of changes of the energy splittings^ 
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FIG. 3: (color online) ESR linewidth (a) and g-factor (b) 
as a function of temperature for the magnetic field applied 
along the three crystallographic axes. The lines in Fig. Ela) 
represent a fit by Eq. 1121 using parameters listed in Table HTl 



Unfortunately, no spectra were shown in Ref. making 
it difBcult to judge where this discrepancy comes from, 
especially, because the spectra were fitted with a single 
Lorentzian line shape in both cases: 
Concerning the uncertainty of the obtained ^-values, 
one has to take into account the strong increase of the 
linewidth with temperature, because the uncertainty of 
the g-value becomes larger as the order of magnitude 
of the linewidth becomes comparable to the resonance 
field of the ESR spectrum (see e.g. Ref. I2l2^ . There- 
fore, we assume the uncertainty in the resonance field as 
5% of the linewidth and obtain the error bars shown in 
Fig. Otb) . Despite these error bars our values at room 
temperature ga = 1.943, gi, = 1.926, and gc = 1.919 dif- 
fer considerably from the ones presented in Ref. |3 (see 
Table0. Here, we would like to emphasize that we inves- 
tigated several samples, which were shown to be of very 
high-quality by clearly revealing the hysteresis at Td in 
the magnetic susceptibility»i^ 

Note that the largest discrepancy for the g-values is 
found for ga, i.e. with the external field applied along 
the a-axis. Interestingly, for this case there is also a 
slight deviation in the temperature dependence of AHa 



in comparison to our data, suggesting that AHa — AHi, 
at room temperature, which we can exclude from our 
data. Hence, the given error for the g-value in Ref. |3 
might have been somewhat underestimated. Moreover, 
none of the calculated sets of g- values obtained by using 
an angular overlap model (AOM) (see Table P) can re- 
produce ga = 2.01 from Ref. ,T&. while the corresponding 
orbital energy levels seem to be in agreement with opti- 
cal data>i^ Instead, the values obtained from the AOM 
in case of isotropic 7r-interaction (model A) describe our 
g-values very nicely. Therefore, we conclude that our g- 
factors correctly reflect the properties of TiOCl and ex- 
clude relevant changes in the crystal-field splitting up to 
room temperature. This is in agreement with direct opti- 
cal measurements of the d-level splittings^^ and the fact 
that x-ray diffraction measurements did not detect signif- 
icant changes of the crystal structure with temperature^ 
With regard to the increase of the g-values towards 
higher temperatures one has to take into account the 
larger uncertainty due to the broadening of the line. In 
principle, however, such a shift could indicate a change of 
the local structure of the Ti04Cl2-octahedra. To decide 
about this possibility, additional structural investigations 
for T > 300 K are desirable. 



IV. g-FACTOR AND CRYSTAL FIELD 
SPLITTINGS OF Ti^+ 



To analyze our experimental ^-factors, we consider 
the local environment of the Ti^^^ ion as a Ti04Cl2- 
octahedron with a strong tetragonal distortion along the 
a-axis. Then we can express the g value parallel (gy) and 
perpendicular {g±) to the direction of tetragonal distor- 
tion as foUowsiSi 



9\\ 



-2-8A||/Ao, g^ = 2-2X^/A' 



(1) 



Here Ay (Aj_) and Ag (A') denote the spin-orbit (SO) 
coupling parameter and the relevant crystal-field split- 
ting, respectively, for the magnetic field applied parallel 
(perpendicular) to the a-axis. Recalling the above discus- 
sion about the increase of uncertainty in the values with 
increasing temperatures (see Fig. OJb)), we will restrict 
the following evaluation to the g- values ga — 1.935(8), 
gb 1.921(5), and g^ = 1.908(5) obtained at 150 K, be- 
cause this temperature is well above Td and the uncer- 
tainty is still quite low. Note, however, that the absence 
of any significant temperature dependence of the factor 



TABLE I: Comparison of the g- values at room temperature 
to experimental results and model calculations of Ref. 
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up to room temperature allows to apply the following re- 
sults in this temperature range with good accuracy. 

Thus, identifying the experimental value of ga with 
g\\ — 1.935(10) and substituting A|| by the isotropic free- 
ion value Ag/fcs = 224 K for Ti^+f2i we derive the en- 
ergy splitting between the ground state and the \d^2_y2) 
level to be Aq = Aj.2_y2 « 2.4(2) eV:^ In compari- 
son to the value A^°'^,^2 = 1-5 eV obtained by optical 

X y 

measurementsiS the value derived from our g-factor is 
too large. Therefore, we have to take into account a co- 
valence reduction of the spin-orbit coupling A|| To esti- 
mate the reduction factor we use the experimental value 
aI?.*^ and obtain Xu/ks ^ 140 K for TiOCl, consider- 
ably smaller than the free-ion value Afi but in agreement 
with literaturetiSi2i This large splitting allows to discard 
the scenario of \dx2-y2) being the first excited state ap- 
proximately 0.34 eV above the ground state (point-charge 
model), in favor of cluster calculations predicting the first 
excited state to be {\dxz) — \dyz)) / ^^^^ 

Additionally, the g-factors in the (6c)-plane can be 
used to estimate the energy of the first exited state in 
TiOCl. In order to simulate the anisotropy of the g- 
factor with the AOM model (see Table QJ, also Kataev 
and coworkers had to consider covalence reduction fac- 
tors, which were chosen to be anisotropic because of a 
presumably stronger covalency of the short Ti-0 bond 
along the a-axis compared with the longer bonds in the 
(6c)-plane. We cannot unambiguously determine the co- 
valence reduction within the (6c)-plane, but we can use 
Aj^ = 140 K and the free ion value Afi = 224 K to obtain 
lower and upper limits for the energy splittings. Starting 
with Aj_ EE A|| = 140 K and = {gi, + gc)/2 « 1.92(1) 
we derive A'^^^^ « 0.3 eV for the energy splitting of the 

doublet \dxz), \dyz) with respect to the ground state. In 
the real structure this doublet splits into the lower an- 
tisymmetric {\dxz) - \dyz))/V2 = |-) (energy Ai) and 
higher symmetric Qd^z) + \dyz))/\/2 = \+) (energy A2) 
state. Using 2/ A' — 1/Ai + I/A2 and the experimen- 
tal value A2 — 0.65(±0.15) eV,^^ we finally obtain the 

lower limit Aj^'*'"'' = 0.2(1) eV. Analogously, we derive 
the upper limit A^^'"' — 0.4(1) eV, narrowing down the 
energy of the first exited state to Ai = 0.2 — 0.4 eV. 
This is in good agreement with the theoretical estimates 
of Ai — 0.25 — 0.3 eV obtained by cluster calculations.^^ 
Thus, by means of ESR we can exclude the degeneracy 
of the first and second excited states in TiOCl, as indi- 
cated by band-structure results,®'^" corroborating the re- 
sults obtained by optics and ARPES measurements >i2iiS 



V. SPIN RELAXATION IN TiOCl 

A. General remarks 

Having identified the character and splitting of the 
ground and low lying excited states via the ^-factors, we 



will now discuss the angular and temperature dependence 
of the linewidth above Td, which provides information 
on the microscopic spin dynamics involving these energy 
levels. The behavior of the ESR linewidth in TiOCl can 
be clearly divided into the three regimes T < 90 K, 90 K 
<T < 250 K, and T > 250 K. In the temperature range 
Tci < T < 90 K the linewidth is almost constant (for 
H\\a,c) or decreases (for H\\b) on increasing tempera- 
ture. This behavior changes at about 90 K together with 
a change of the linewidth anisotropy (see Fig. [7^ to the 
monotonous increase for all directions. At high temper- 
atures T > 250 K a strong additional increase of AH{T) 
dominates this saturation behavior (see Fig. 0) and the 
anisotropy of the line vanishes (Fig. EJ. We attribute 
these different regimes to the competition of relaxation 
mechanisms prevailing at different temperatures, which 
will be discussed in detail in the following. At first we 
have to single out the relevant interactions which drive 
the relaxation in TiOCl: 

Single-ion anisotropy is absent for Ti'^+ (S'=l/2). 
Other sources of line broadening such as dipolar inter- 
action or hyperfine coupling are negligible as a result 
of the large isotropic exchange J/fce = 660 K. Taking 
into account the average distance between the Ti-ionai 
of about 3.355 A and the value of the hyperfine con- 
stant Ar^i3+ = 6 • 10~^ cm~^ (see Ref. [25f) we can esti- 
mate the contribution to the linewidth from these sources 
as 10^^ Oe and 10~^ Oe, respectively. The minor im- 
portance of these interactions for the linewidth broaden- 
ing in low-dimensional systems with a strong exchange 
coupling has been discussed in detail e.g. by Pilawa et 
al.'^^ for CuCeOs and Yamada et al.^'' for NaVsOg. A 
larger contribution to the ESR linewidth could be ex- 
pected for the anisotropic Zeeman interaction in case of 
different Ti'^+ sites in adjacent layers»2S, However, this 
broadening strongly depends on the value of the reso- 
nance field i?res- At X-band frequency used in our ex- 
periment (Hros ~ 3 kOe) the resulting contribution is 
less than 1 Oe for any reasonable choice of parameters 
(e.g. an interlayer couphng Jintcr ~ 0.05 • J fRef. [lo|) and 
Ag ~ 0.3). 

The remaining relevant contributions stem from 
the anisotropic exchange interactions. Conventional 
estimations^® of their magnitude result in values at least 
two order of magnitudes higher as for the other sources of 
line broadening. Note, however, that the applicability of 
such estimations for low-dimensional systems like TiOCl 
has been questioned, recentlyi22i22i2i Hence, one has to 
analyze carefully each spin system under consideration 
on a microscopic level. 

In general, anisotropic exchange interactions arise due 
to virtual hopping processes of electrons or holes between 
two interacting ions via the bridging diamagnetic ions in 
combination with the SO coupling Hso = XLaSa- Here 
La and Sa denote the orbital and spin momentum of 
the magnetic ion, respectively. For a more detailed dis- 
cussion of the exchange interactions we refer to Refs. 3^ 
and I33I Now, we will discuss the symmetric part of the 
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anisotropic exchange (AE) interaction as well as the anti- 
symmetric Dzyaloshinsky-Moriya (DM) exchange inter- 
action and their importance for the spin relaxation in 
TiOCl. Their influence on the spin relaxation will be an- 
alyzed in terms of the so called moment method: In the 
case of sufficiently strong exchange interaction the ESR 
linewidth can be obtained from the second moment 



Mo 



([5+,5-]) 



of the ESR line via the relatio: 



,34 



AH 



1 M2 

5MB J ' 



(2) 



(3) 



where /ie denotes the Bohr magneton. This method 
is well established only in the high-temperature limit 
T > J/fce. In this limit the second moment is tem- 
perature independent and can be expressed via the pa- 
rameters of the microscopic spin Hamiltonian Ti-mt^ 
As the exchange constant is quite large in TiOCl, we 
do not reach the high-temperature limit in our exper- 
iment. However, the anisotropy of the linewidth yields 
the relative values of the exchange parameters even at low 
temperatures:'^'' In the following we determine the lead- 
ing components of the symmetric AE tensor and of the 
DM vector. Inserting these parameters into Eq. Q yields 
the angular dependence of the linewidth which then can 
be compared to the experimental observations. 



far the most relevant ones. The first exchange process 
between neighboring Ti ions in the chain (depicted at 
the left side of the upper panel) with an electron trans- 
fer between \dx2_y2) and |c?2,y)-orbitals becomes impor- 
tant as a result of the strong cr-bounding between the 
titanium d- and the oxygen p-orbitals. The importance 
of the second intra-chain AE process (at the right side 
of the upper panel of Fig. ^ via |— )-orbitals is due to 
the small energy Ai of the involved exited state. Note 
that only the exchange paths between the exited |— )- 
states are shown in the second case, since the exchange 
path between the ground |(ia;y)-orbitals has already been 
discussed elsewhere]^ In the lower panel of Fig. 0] the 
dominating exchange paths of the inter-layer AE are pre- 
sented: the left one being between |— ) and |da;i/)-orbitals 
and the second between |-|-) and |(ia:j,)-orbitals. These 
processes are of the same order of magnitude as the intra- 
chain exchange and cannot be neglected. 

The non-zero elements of the exchange tensors can 
be determined via the SO operators included in this 
processi^SiSiSi For the first intra-chain AE process (via 
the |fi3.2_j^2)-orbital) the exited state is connected to 
the ground state dxy of the same Ti ion via SO cou- 
pling with only one nonzero matrix element, namely 
{dx.2_y2\Lz\dxy) = —2i. Following Ref. Is^ all AE pro- 
cesses via this level contribute to Jzz only. Taking now 
into account the relation for the diagonal components of 
the AE tensor ^ Jaa = we obtain for this process 



B. Symmetric anisotropic exchange 

The importance of symmetric anisotropic exchange for 
low-dimensional systems has been emphasized in a recent 
study by Eremin and coworkers for Q;'-NaV205i^ It was 
shown that it is necessary to take into account both the 
exchange geometry and the energies of the involved or- 
bital states in order to obtain reliable estimates of the AE 
on a microscopic basis. All electron transfers processes 
both between the exited states on two sites and between 
the exited and ground states of the involved ions must 
be considered. Then, contributions of possible pathways 
to the AE can be calculated using the respective values 
of hopping integrals and energy splittingsn^ Such micro- 
scopic estimations are beyond the scope of this paper, 
but we will discuss qualitatively the relevant exchange 
processes in TiOCl based on the analysis of Ref. [s^ 

Crystal-field field splittings of the relevant exited states 
have been already estimated above (see Sec. IIV|) and in 
order to illustrate the exchange geometry via these or- 
bital levels we show the charge-distribution pictures for 
the pathways of AE for the intra-chain and inter-layer 
exchange in Figs.01a) and^b), respectively. The inter- 
chain AE within one layer is not effective here because of 
the orthogonality of the ground-state orbital with respect 
to the direction of the exchange. 

In analogy to the estimations made in Ref. 36 we ar- 
gue that the pathways of AE shown in Fig. 0] are by 



_2 7(l^ -r» 



^Jyy 



(4) 



All other AE processes which make a considerable con- 
tribution to the linewidth in TiOCl involve the |— ) or |-|-)- 
orbitals (|±) = {Id^z) ± Myz)}/%/2), which are connected 
to the ground state orbital \dx-^ via the matrix elements 



{dxz\1^x\dxy) 



i and {dyz\Ly\dxy) 



The resulting 



nonzero elements Jxx and Jyy of the AE tensor have the 
same magnitude because of symmetry reasons and they 
have the same sign, because the expression for Ja/s (see 
Ref. depends on the square of the orbital momentum. 
Thus, we can write the AE tensor for these processes as 



= -2Jir'''^'^ = -2J^Jy'''^^^>. It becomes clear 
that the maximal component of the anisotropic exchange 
tensor is Jzz- Therefore, we would expect the maximal 
linewidth for H\\z in agreement with the experimental 
data for T > 90 K. 

Consequently, we can describe the resulting angular 
dependence of AH in terms of the moment method in 
analogy to CuGeOg (Ref.lH) 



r(|->,l+» _ 



t{|->.|+)) 



AHae{T, 9) = Kae(T){1 + cos^ 



(5) 



where 9 is the polar angle of H with respect to the a = —z 
axis and Kae (T) is proportional to the strength of the 
AE interaction. This parametrization does not allow to 
obtain the exact values of the anisotropic exchange pa- 
rameters, but it is valid for all temperatures and, hence, 
we will apply it to describe our data using Kae as a fit 
parameter. 
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FIG. 4: (color online) Schematic pathway of the symmetric 
anisotropic exchange between Ti ions in TiOCl. Big spheres 
denote Ti ions, small spheres - O and CI ions. Upper frame: 
The most relevant intra-chain (|| b-axis) exchange paths. Left 
one - between the ground state \dxy)-orhital on one site and 
the exited |ci2.2_j^2)-orbital on the other site, second one - 
between two exited | — )-orbitals on both sites. Lower frame: 
two dominating exchange paths of the inter-chain AE between 
a Ti ion in the \dxy)-sta.te and the nearest-neighbor Ti ion 
from the adjacent chain in the exited states | — ) (left) or |+) 
(right) . 



Concerning the temperature dependence of the ESR 
linewidth produced by AE exchange interaction, clear 
theoretical predictions do only exist in two limiting 
cases: (i) For the high-temperature regime T > J, 
the linewidth approaches the result of the Kubo-Tomita 
theory34 Ah'xe^\t) Kae{oo) cx J^^/J = const, 

and (ii) the result AiJ^^<<"'^ (T) cx {Kae{oo)/J) ■ T in 
the case T <^ J for the S — 1/2 quantum antifer- 
romagnetic chain^ To model the crossover regime we 



will use the following empirical expression which has al- 
ready provided a successful description for several low- 
dimensional systems like CuGeOsj^ LiCuV04f?^ and 
Nai/gVaOs (Ref.lH: 

AHae{T) = Kae{T) = Kae{^) • e-TC. (6) 

Thus, the fit parameters to describe the contribution of 
the AE are Kae{oo),Ci,C2- 

C. Dzyaloshinsky-Moriya interaction 

The contribution of the DM interaction to the ESR 
line broadening in one-dimensional systems is a heavily 
debated topic at the moment i2L^2i2ii^ Using the well- 
established high-temperature Kubo-Tomita approach^M 
the DM interaction was considered to be the dominat- 
ing relaxation mechanism in spin-chain compounds like 
e.g. NaV2 05i^ However, the applicability of this ap- 
proach was questioned in a field theoretical one by Os- 
hikawa and Affleckf^ arguing that the contribution of the 
DM interactions is strongly overestimated by the Kubo- 
Tomito approach. A recent conclusive experimental con- 
formation of the results of Oshikawa and Affleck was re- 
ported by Zvyagin et al. , who unambiguously showed 
that the importance of the DM interaction in antiferro- 
magnetic spin- 1/2 chains has to be considered with spe- 
cial care4i While the DM interaction was not explicitly 
discussed by Kataev and coworkers for TiOClr^ Kato et 
al. considered it to be the main source of line broaden- 
ing in the isostructural system TiOBr by referring to the 
Kubo-Tomita approach4^ 

In general, the structure of TiOCl allows for a contri- 
bution of the antisymmetric DM interaction 

HI;^, - d,,, . [S,; X S,] (7) 

between the Ti spins Si and Sj via an intermediate dia- 
magnetic ioni22i^ The DM vector is an axial vector 
perpendicular to the plane spanned by the Ti spins and 
the ligand ion, determined by Aij = d^'^ • [xiu x rij;], 
where the unit vectors wa and rij-; connect the spins i 
and j with the bridging ion /, respectivelyiii When the 
point bisecting the straight line connected two interact- 
ing ions is not a center of inversion, which is the case for 
TiOCl or TiOBr^ one can expect that dy ^ Oi2£ Kato et 
al.'^'^ argued that the maximal component of the DM vec- 
tor in TiOBr is parallel to the a axis, accounting for the 
maximal linewidth along the a axis in the temperature 
regime T^i < T < Tc2. Note, that in TiOCl the max- 
imum linewidth in this temperature regime is observed 
for II 6 (Fig. Hand Ref.[li). 

In what follows, we will analyze the DM interaction 
in TiOCl considering every next-neighbor bond of the Ti 
ion on basis of the structural data of Shaz et al. at room 
temperature.^ All bonds of the Ti ion together with the 
corresponding DM vectors are shown in Fig. [S] Only in- 
teractions of Ti ions in the same layer give rise to the an- 
tisymmetric exchange in TiOCl because of the existence 
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df')(-1A 



("^^(1,0,0) 



dJ0,1,0) 




r(o) 



FIG. 5; (color online) Next-neighbor bonds of the Ti ion to- 
gether with the associated parameters of the antisymmet- 
ric anisotropic exchange. dt, di'^' and di'^^' denote the 
Dzyaloshinsky-Moriya parameters for the exchange along the 
a-axis, along the fe-axis via the oxygen ion and along the b- 
ajcis via the chlorine ion, respectively. Components of the DM 
vectors are given in the depicted cartesian coordinate system 
{b,c,a}. Big spheres denote Ti ions, small spheres - O and CI 
ions. 



of an inversion center between the Ti ions from adjacent 
layers. The two remaining contributions arise from the 
chains of the Ti ions along the b and a directions. The 
first one, which results in a component of dij in the a 
direction, has been considered in Ref. |^ as the domi- 
nating source of the line broadening. However, we would 
like to point out that there are two different bridging 
ions (Cl~ and 0^~) leading to DM vectors with opposite 
sign in this case. Although the two paths are asymmet- 
ric (nxi-ci ^ 2.393 A, nxi-o ^ 2.187 A at T = 295 K)i 
and, hence, lack inversion symmetry, one can assume that 
the opposite DM vectors will partially compensate each 
other. If we denote the respective DM parameters 
as d^a^^ and c?i^'^ for the exchange via the and Cl~ 
ions, respectively, only its difference Ada = d'p^ — di^'^ 
will give rise to the ESR line broadening and can be de- 
tected experimentally (see e.g. the discussion about the 
cancelation of DM interaction in LiCuV04)^^. Looking 
now at the contribution of the inter-chain DM interaction 
(between two neighboring Ti^+ sites along the a axis via 
the oxygen ion lying in the same (ac)-plane) , we can con- 
clude that the corresponding DM vector is pointed along 
the 6-axis (see Fig.EIl. 

The general expression for the M-? due to the DM 
interaction has been given in Ref. In the case of 
TiOCl only two intra-chain and two inter-chain contri- 
bution must be taken into account by the calculation 
of the respective ESR line broadening. Estimation of 
the "geometrical factors" G^'-* — [n.a x nji] yields: (i) 



±0.98 « 1 and Gr ^ = ±0.99 w 1 for the 
intra-chain exchange via the O^^ and Cl~ ions, respec- 
tively. Therefore, we obtain the following expression for 
the second moment of the DM interaction in the crystal- 
lographic system: 



2 {e,^)(x (4-i/2)(i + 

+ (Ada-l)(l 



sin ^cos^iyj)- 



(8) 



where 9 and if are the polar and azimuthal angles of 
H with respect to the a axis. Finally, the ratios of the 
linewidth along the three crystallographic axes read 



AHh : AHn : AHr 



M2(f,0) :M2(0,^) :M2(|,f) 



l+2{2 Ada /djy 
l + (2Ada/di,Y 



(9) 



Simplifying this expression for the case Ada/db — > 0, one 
gets AHb : AHa : AHc = 2:1:1 and the angular 
dependence 



AH 



(Ada/db^O) 
DM 



(r,0,(^) =i^Z5M(r)(l+sin20cosV), (10) 



where if is the azimuthal angle of H in (6c)-plane with 
respect to the b axis and Kdm(T) is proportional to the 
strength of the DM interaction. 

Regarding the temperature dependence of AHjjm (or 
Kdm{T)), we will use the result obtained by Oshikawa 
and Affleck AH^j^^'^\t) cx J'^/T'^ obtained for the case 

of a staggered DM interaction H^^^lj = J2i ' I^* ^ ^j+i] 
with dj = (-l)*d for T < JM Assuming that this tem- 
perature dependence also holds for a uniform DM inter- 
action along the chain as in our case, we will apply the 
power law 



AHdm{T) = Kdm{T) = Kdm{^) { ^ 



(11) 



for the inter-chain exchange G 



(O) 



±0.501 w 1/2, (ii) 



to fit the experimental data using K]j]\i{oo) as a fit pa- 
rameter. An analytical expression for the crossover be- 
havior from this power-law (valid for T <C J = 660 K) to 
the constant high-temperature value of the Kubo-Tomita 
approach has not been derived up to now. Therefore, 
we extrapolate the power law up to T = J and iden- 
tify Kdm{T — J) = Kdm{oo) in order to compare the 
experimental values to the theoretical estimates of the 
Kubo-Tomita approach. We would like to recall that the 
ratio of AH along different axes in Eq. ^ evaluated in 
the high-temperature limit does not depend on the form 
of the temperature dependence. 



VI. ANALYSIS AND DISCUSSION OF THE ESR 
LINEWIDTH 

Starting out with the superposition of the angular and 
temperature dependence for the AE and DM interac- 
tions, it is possible to describe the anisotropy and the 
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temperature dependence of the linewidth for Td < T < 
250 K, but the change of curvature (Fig.EJ) for T > 250 K 
clearly shows that an additional relaxation channel dom- 
inates at higher temperatures. Before we discuss the low- 
temperature data in detail, we shortly comment on pos- 
sible reasons for this high-temperature behavior. 

In order to determine the temperature dependence 
more accurately up to 500 K, we additionally performed 
measurements of crushed single crystals (see Fig. El . This 
was necessary because of the fact that the single crystals 
of TiOCl are thin platelets of small mass and that the 
linewidth above room temperature is already very large. 

It turns out that the strong increase of the ESR 
linewidth with temperature can be very well accounted 
for by adding an exponential term ifexp • e"^^^"/*^^"^ to 
the temperature dependence of the anisotropic exchange 
interactions: 



AiJ(T) = Kdm{^) ■ {^y 



+Kae{oo) ■ e 
e 



' T + C2 



(12) 



7^ . -AESB/kBT 

^ ^ cxp ' 



The resulting fit is shown as a solid line in Fig. El yielding 
^ESR = 0-28 eV. The exponential nature of the addi- 
tional increase is highlighted in the inset of Fig. El where 
the reduced linewidth data Ai?rcd are plotted after sub- 
traction of the contributions of the AE and DM interac- 
tions (dashed hne in Fig.EJ. 

An additional relaxation channel via thermally ac- 
tivated charge carriers might cause an exponential in- 
crease of AH{T), as it has been discussed for doped 
manganites^ and at the metal-to-insulator transition 
in /3-Nai/3V205i^ In both cases the leading contribu- 
tion to the temperature dependence is determined by the 
Arrhenius law of the conductivity . The cor- 

responding temperature behavior of the dc-conductivity 
which could be obtained from dielectric measurements4& 
In an Arrhenius representation one can extract an acti- 
vation energy Aa/2 « 0.31 eV (see inset of Fig. EJ."'^ 
Although this value is similar to the one obtained from 
the ESR linewidth, it is by far too small compared with 
the experimental gap value of about 2 eV observed by 
optical spectroscopy.^'* Therefore, such a scenario ap- 
pears rather unlikely. Alternatively, the exponential in- 
crease can be interpreted analogously to the case of the 
one-dimensional magnet CuSb2 06. In this compound a 
similar temperature dependence of the linewidth with 



A 



(CuSb^Oe 
ESR 



0.13 eV was observed and explained in 



terms of a thermally activated dynamic Jahn- Teller (JT) 
process."*^ However, a rigorous theoretical treatment of 
this effect has not yet been undertaken. Moreover, we 
would like to mention that the obtained value Aesr ~ 
0.28 eV is very close to the energy 0.3(1) eV of the first 
excited state. This might indicate the involvement of ex- 
ited orbital states in this relaxation process. To finally 
decide about the origin of the high-temperature relax- 
ation, however, detailed structural studies in this tem- 
perature region are necessary. 




FIG. 6: (color online) Temperature dependence of the ESR 
linewidth for a powder sample together with a fit by Eq. 1121 . 
The obtained fit parameters are given in Tablejni Aii'ox(r) = 
AHae{T) + AHdm{T) is shown as dashed line. Inset: re- 
duced contribution A//rcd(r) = AH{T)-AHa^{T) compared 
to conductivity as Arrhenius plot. 



Fixing the value Aesr — 0-28 eV, we now proceed to 
describe the anisotropic temperature dependence of the 
linewidth for the main orientations of the single crystal. 
Note that Aesr should not depend on the orientation 
of the magnetic field with respect to the crystal axes, 
justifying the further use of this value for the single crys- 
tal. The resulting fit curves are shown in Fig. El and 
the obtained fit parameters are given in Table ^ The 
agreement between fit and data below 250 K is excellent 
(see also Fig. [71), but at higher temperatures deviations 
are clearly visible for H\\ a and H\\b. The anisotropy in- 
ferred from the AE and the DM interaction below 250 K 
is somewhat larger than the observed one. The gradual 
suppression of the anisotropy with increasing tempera- 
ture may result from the thermal occupation of higher 
lying c?-levels, which is also in agreement with the dis- 
appearance of the anisotropy of the g-tensor (Fig. EJ. 
A similar effect was observed at the transition from a 
cooperative static JT-effect to a dynamic JT-phase in 
(La:Sr)Mn03-2^ 

Let us turn to the discussion of the anisotropic ex- 
change contributions which dominate the relaxation be- 
low 250 K. Using the obtained fit parameters we addi- 
tionally fitted the angular dependence of the linewidth 
data for the single crystal in the crystallographic (ab)- 
plane (Fig. O by using 



Aff(°^)(r,( 



_ A _ ^ , 

e 



+ Sae ^^AE 



sm 

_c: 

T+C; 



(1 



COS 
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FIG. 7: (color online) Temperature dependence of the ESR 
linewidth in TiOCl for the magnetic field applied parellel to 
the three crystallographic axes. The fit curves were obtained 
by the sum of AHae{T) and AHdm{T) (Eqs. ©.CJ) with 
the parameters listed in Table|n| The dashed and dotted lines 
represent the contributions from the AE and DM interactions, 
respectively. 



500 



400 



O 300 



X 

< 



200 



100 





Hllb 


Hlla 

* 


— 








TiOCI 








/ A 150K ^ 
/ o 100K 






^ 70 K 






^nOOOOoOn 

^^^.^ 












^^^i^A A A A A A A A A^i-^f*^ 





60 


120 180 



angle (deg) 

FIG. 8: (color online) Angular dependence of the ESR 
linewidth AH in TiOCI for the magnetic field applied within 
the (a6)-plane at different temperatures between 70 and 
150 K. The solid lines represent the fit by the Eq. 1131 . where 
the fitting parameters are taken from Table ITTl 



(13) 

where 9 denotes the angle in the (a5)-plane with re- 
spect to the a axis. Here, we took into account only 
the DM contribution along the crystallographic 6-axis 
(see Sec. IV Cll . Concerning the AE interaction, we 
had to introduce the additional fit parameter Sae ~ 
k'xJ"'\oo)/[2 ■ k'xJ''\oo)] = 1.05 which indicates the 
deviation from the theoretically expected ratio of 1, if 
only the AE paths described above are taken into ac- 
count (i.e. Jzz — —'i.Jxx = ~'^Jyy)- The fact that 
Sae = 1-05 can be explained by small contributions of 
the other relaxation processes (see Sec. IV At . Thus, we 
were able to corroborate the validity of the fit parame- 
ters given in Table ^1 by a consistent description of the 
temperature and angular dependence of the linewidth. 
Moreover, we find a good agreement of the ratios of the 
obtained high-temperature fit parameters for the DM 
interactions ^^^^^^(oo) : if^^|"^(cx)) : k\^j^I^'> (oo) = 
ATjf ^ : AH^^'^ : AH^^'^ = 1.72 : 1.04 : 1 with the 



theoretically expected ratio 2:1:1. 

Looking at the corresponding contributions of AE and 
the DM interactions shown in Fig. [T] it becomes clear 
that the dominant relaxation mechanism for T > 90 K 
is the AE, while the DM interaction takes over for Td < 
T < 90 K. This competition is nicely evidenced by the 
corresponding orientation dependences and the crossover 
at about 90 K (Fig. (S)). However, significant contribu- 
tions of the DM interactions can already be anticipated 
below 135 K where the linewidth for H \\ b already be- 
comes larger than the one for H \\ c. Here, we have 
to emphasize that our analysis of the DM interaction is 
based on the room-temperature structure and does not 
take into account a possible structural phase transition 
at Tc2- Since the linewidth data does not reveal a dis- 
continuity at Tc2 but a smooth crossover, we conclude 
that the structural changes do not significantly alter the 
involved relaxation processes. 



VII. SUMMARY 



TABLE II: Parameters determined from fits on the temper- 
ature (Figs. 13161711 and angular (Fig. (HJ dependencies of the 
linewidth by the Eq. 1121 and 11311 . respectively. The pa- 
rameters Ci = 129.12 K, C2 = -38.1 K, Aesr = 0.28 eV, 
Js'exp = 0.87 MOe are assumed to be isotropic. Parameter 
Sae used by the fit of the angular dependence of AH in the 
(o6)-plane is equal to 1.05. 





/Cab (00) [Oe] 


Kdm(oo) [Oe] 


H\\ a (single crystal) 


1429 


1.397 


H\\b (single crystal) 


765 


2.319 


H\\c (single crystal) 


930 


1.344 


crushed single crystal 


990 


1.344 



In summary, we have investigated the temperature de- 
pendence (T < 500 K) of the ESR linewidth AH and the 
g value in TiOCI. From the g values we derive the energy 
of the first excited state as Ai = 0.2 — 0.4 eV, in good 
agreement with theoretical estimations. Furthermore, we 
describe the angular and temperature dependence of the 
linewidth as a competition of the anisotropic exchange 
interactions and an additional exponential increase for 
T > 250 K higher temperature that might be related to 
thermally activated lattice fluctuations. We could show 
that the line broadening is dominated by the symmetric 
anisotropic exchange for 90 K < T < 250 K which pro- 
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duces the maximal linewidth along the a direction, while 
the antisymmetric DM interaction leads to the crossover 
at about 90 K with the maximal linewidth along the h 
direction. 
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